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Abstract 

The degree of synchronization and the amount of dynamical cluster formation in 
electroencephalographic (EEG) signals are characterized by employing two order 
parameters introduced in the context of coupled chaotic systems subject to exter- 
nal noise. These parameters are calculated in EEG signals from a group of healthy 
subjects and a group of epileptic patients, including a patient experiencing an epilep- 
tic crisis. The evolution of these parameters shows the occurrence of intermittent 
synchronization and clustering in the brain activity during an epileptic crisis. Sig- 
nificantly, the existence of an instantaneous maximum of synchronization previous 
to the onset of a crisis is revealed by this procedure. The mean values of the or- 
der parameters and their standard deviations are compared between both groups of 
individuals. 
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1 Introduction 



In recent years there has been substantial interest in the description of neural 
processes and electroencephalographic (EEG) signals in the context of nonlin- 
ear dynamics and deterministic chaos [1,2,3]. Nonlinear dynamics opens new 
windows for the understanding of EEG signals since the neural activity in- 
volves nonlinear mechanisms at the microscopic level. The complexity of the 
behavior of neural systems reflects the effects of those underlying nonlinear 
mechanisms. In particular, the occurrence of synchronization in EEG signals 
associated to various brain disorders, such as epilepsy and migraine, has been 
the focus of much attention [4,5]. Synchronization phenomena are expected to 
play a major role in establishing the communication between different regions 
of the brain [6,7]. 

The phenomenon of chaotic synchronization takes place when the phase space 
trajectories of two or more coupled chaotic systems, initially evolving on dif- 
ferent attractors, eventually converge to a common trajectory [8,9]. A related 
phenomenon is dynamical clustering; it consists of the formation of differenti- 
ated subsets of synchronized elements in a coupled system [10]. Experimental 
observation of dynamical clustering in coupled electrochemical systems has 
been reported [11]. 

This article presents a procedure that allows to characterize the amount of 
synchronization and clustering occurring on coupled chaotic oscillators sub- 
ject to common noise and applies these concepts to EEG signals from healthy 
subjects and epileptic patients. In Section 2, a theoretical model consisting of 
a set of globally coupled chaotic Rossler oscillators and subject to an external 
noise is introduced. Synchronization and the formation of clusters, i.e., syn- 
chronized domains, under the effect of noise are studied in this model. Two 
order parameters are introduced in order to characterize the synchronization 
and the formation of clusters, i.e., synchronized domains, under the effect of 
noise in this system. We calculate these order parameters for different values 
of the coupling strength between the oscillators. The data base of EEG signals 
used in this study is described in Section 3; it consists of a group often healthy 
subjects and a group of ten epileptic patients, including one EEG signal of a 
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patient experiencing an epileptic crisis. In Section 4, we perform the analysis 
of synchronization and cluster formation in the EEG signals from the data 
base. The analysis yields quantitative indices that characterize the changes of 
complexity in the brain activity of healthy and epileptic subjects. Conclusions 
are presented in Section 5. 

2 Globally coupled Rossler oscillators subject to external noise 

The Rossler equations constitute a paradigmatic dynamical model that ex- 
hibits chaos. For some parameter values, these equations give rise to a chaotic 
attractor possessing fractal properties. In order to introduce the concepts used 
in the characterization of EGG signals, we consider a theoretical model consist- 
ing of system of globally coupled Rossler oscillators and subject to an external 
noise described by the following equations [12]: 

ii = -Vi - Zi + e(x - Xi) + iiit) 

y< = Xi- ayi + e(y - yi) (!) 
Zi = b - czi + XiZi + e(z - Zi) , 

where Xi, yi and Zi are the state variables of oscillator i; with % — 1, . . . , N; N 
is the number of oscillators; e represents the strength of the coupling amount 
the oscillators; is the source of noise acting on oscillator i, with mean 
value (£i(t)) = and correlation = 2S5(t — t')8ij\ and S is the 

amplitude of the noise. The mean value of the variables Xi is given by 

1 N 

z(0 = T7l>W, (2) 
iv i=l 

with similar expressions for y and z. The addition of noise in Eq. (1) is relevant 
because we wish to establish a valid comparison between the behavior of this 
dynamical model and that of physiological systems which usually present noise 
and fluctuations, as in the case of EEG signals. The values of the parameters 
a, b and c employed here are a = b = 0.2, and c = 4.5 for which the individual 
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Rossler oscillators are chaotic. 

A cluster is defined as a subset of the population of oscillators which are 
synchronized among themselves. However, the presence of noise in a system 
of interacting elements usually prevents perfect synchronization. In practice, 
we consider that a pair of elements % and j belong to a cluster if the distance 
dij in phase space between them, defined as 

dij = \j{x t - xj) 2 + (yi - yj) 2 + 0; - zj) 2 , (3) 



is less than a threshold value 7, i.e., d^ < 7. The choice of 7 should be 
appropriate for achieving differentiation between closely evolving clusters. 

To characterize the amount of synchronization and the formation of clusters 
in a system of coupled chaotic elements such as Eq.(l), two different order 
parameters can be introduced [13,14]. The first one, r(t), is defined as the 
fraction of pairs of elements which are separated by a distance d^ < 7 
at time t, 

1 N N 

r(t) = N{N _ 1) j:T i eh-d ij (t)}, (4) 

where 9(x) is the step function such that 6(x) = for x < , and 6(x) = 1 
for x > 0. The second order parameter, is the fraction of elements i which 
at time t have at least one other element j located at a distance d^ < 7, 

1 N N 

«W = l-jyE II (5) 

i=l j=l,j^i 

The last term in Eq.(5) is the fraction of elements whose distance to any 
other element is greater than 7. Total synchronization corresponds to r(t) = 
s(t) = 1. On the other hand, if all the population of oscillators is distributed in 
clusters, we have s(t) = 1 while r(t) < 1, since the clusters may be separated 
[12]. If some elements are in clusters but others remain loose, then r(t) < 
s(t) < 1. 
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Fig. 1. Order parameters r(t) and s(t) as function of coupling e, at time t = 5000; 
with N = 1000. 

System (1) with N = 1000 and noise of amplitude S = 10~ 2 was integrated 
by using a fourth-order Runge-Kutta method with time increments At = 
10 -2 , starting from random initial conditions. Noise was introduced through 
a reliable random number generator. 

Figure 1 shows the order parameters r(t) and s(i) as functions of e at time 
t = 5000. The threshold for cluster differentiation was set at 7 = 10~ 3 , a value 
for which both r(t) and s(t) reach statistically stationary values. For small 
values of e, both r(t) and s(t) are less than 1, indicating that not all oscillators 
are forming clusters. In the interval e G (0.007,0.090), s{t) increases more 
rapidly than r(t); this indicates that the oscillators tend to synchronize and 
to form clusters. For e = 0.091, it is found r(t) = s(t) = 1, and therefore the 
elements are totally synchronized and forming a single cluster. For e > 0.133 
both order parameters definitively reach their maximum values r(t) = s(t) = 1 
and the system is totally synchronized. 

3 Data base of EEG signals 

The EEG data base used in this study consist of records from 10 healthy 
subjects and 10 epileptic patients. One of the epileptic records belongs to an 
individual who experienced a spontaneous convulsive crisis during the EEG 
recording. The record of the EEG signal from each individual was done over 18 
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Fig. 2. (a) Channel C4 of the EEG signal from healthy subject, (b) Channel C4 of 
the EEG signal from an epileptic patient without crisis, (c) Channel C4 of the EEG 
signal from an epileptic patient experiencing a crisis. The flat extremes correspond 
to saturation of the electroencephalograph. 

channels connected to scalp electrodes according to the international system 
10 — 20. The potentials were measured with respect to a reference level con- 
sisting of both ears short-circuited. The signal was digitalized at a sampling 
frequency of 256 Hz and A/D conversion of 12 bits, and filtered to bandwidth 
between 0.5 and 30 Hz. 

Figure 2 shows the records from the channel C4 of the EEG signals correspond- 
ing to a healthy subject, an epileptic patient without crisis, and an epileptic 
patient just before and during a crisis, respectively. 

4 Results 

In order to characterize the synchronization and cluster formation in the EEG 
signals for each individual, we have employed the order parameters r(t) and 
s(t) introduced in Sec. 2. At time t (ms), the values of r{t) and s(t) were 
calculated for each EEG signal comprised of 18 channels. The value of 7 used 
corresponded to 40% of the standard deviation of each EEG signal. This value 
of 7 is enough to discriminate different clusters in the EEG signals. 

Figure 3 shows the order parameters r(t) and s(t) for a healthy individual. The 
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Fig. 3. Order parameters r(t) and s(t) for the EEG signal from a healthy subject. 

parameter r(t) exhibits relatively small values indicating desynchronization 
among the different channels. The parameters (t) reaches higher values that 
suggest the presence of some clusters. The maximum value reached by s(t) is 
0.78. The analysis of the EEG signals from the others healthy subjects shows 
similar results to those shown in Figure 3. 

Figure 4 shows the order parameters r(t) and s(t) for a epileptic patient ex- 
periencing a spontaneous epileptic crisis. The corresponding EEG signal is 
shown in Figure 2(c). The intervals where the amplifier was saturated were 
excluded for the calculation of r(t) and s(t). Note that the parameter r(t) 
reaches a maximum value r(t) = 0.51 just two seconds before the beginning 
of the epileptic crisis, while at the onset of the crisis this parameter has a 
value of r(t) = 0.39. Several maxima of r(t) are observed during the crisis. 
These maximum values of r(t) are greater than those reached in EEG signals 
from healthy subjects. Similarly, the parameter s(t) in Figure 4 also exhibits 
a maximum value of 0.96 at two seconds before the onset of the crisis. In ad- 
dition to this precursory maximum, the parameter s(t) shows several maxima 
during the epileptic crisis. Most of these maxima correspond exactly to the 
maxima of the parameter r(t). This can be interpreted as the occurrence of 
intermittent synchronization in the brain activity during an epileptic crisis. 

Figure 4 reveals the occurrence of a synchronization maximum at two seconds 
prior to the beginning of the epileptic crisis. However, this is not evident from 
the observation of the EEG signal from the patient, shown in Figure 2.c. This 
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Fig. 4. Order parameters r(t) ands(i) for the EEG signal from the epileptic patient 
experiencing a crisis. The dotted vertical line indicates the onset of the crisis. 

feature shows the potential that the use of nonlinear dynamics techniques has 
in the analysis of physiological signals and how these techniques can contribute 
to the prediction and diagnostics of some pathologies. 

The parameters r(t) and s(t) calculated from the EEG signals from healthy 
subjects and epileptic patients not experiencing crisis show a similar behavior. 
This might be due to the fact that these patients were being treated with drugs. 

Table 1 

Mean values and standard deviations of parameters r(t) ands(t). 
Healthy Subjects Epileptic patients Epileptic patient 
without crisis with crisis 



r 0.075 

Ar 4.12xl0" 4 

s 0.607 

As 4.28xl0~ 3 



0.076 

2.62xl0" 4 
0.604 

3.98xl0~ 3 



0.110 

4.06xl0~ 3 
0.646 
6.14xl0~ 3 



Table 1 displays the mean values r, s, and standard deviations Ar and As of 
the parameters r(t) and s(t) over time and over the population of healthy sub- 
jects, epileptic patients without crisis at the time of recording and the epileptic 
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patient experiencing a crisis. In Table 1, note that the mean value r clearly 
differentiate an epileptic patient with crisis from the other groups. The value 
of r does not unambiguously distinguish between epileptic patients without 
crisis and healthy subjects; however the value of the standard deviation Ar in 
the former group is about twice the value corresponding to healthy subjects. 

5 Conclusions 

Globally coupled systems, such as Eqs. (1), constitute theoretical models con- 
taining the essential ingredients that allow the emergence of synchronization 
and the formation of clusters in dynamical systems consisting of many inter- 
active elements. The central hypothesis of this article is that the brain can 
also be described as a dynamical system formed by mutually coupled nonlin- 
ear elements (neurons), and therefore the collective propierties of the brain 
are susceptible to be analyzed through nonlinear dynamics techniques. With 
this aim, we have quantified the amount of synchronization and clustering in 
EEG signals by using order parameters previously introduced in the context 
of coupled chaotic oscillators. 

In the case of the epileptic patient experiencing a crisis, the order parameters 
r(t) and s(t) show a series of instantaneous maxima in time which are absent 
in the evolution of those parameters in both healthy subjects and epileptic 
patients without crisis. The maxima of r(t) indicate instantaneous increases 
in synchronization. During an epileptic crisis the mean value r is greater than 
those values corresponding to healthy subjects and to epileptic patients with- 
out crisis. The standard deviation Ar is an order of magnitude less in the 
epileptic patient in crisis with respect to healthy subjects and to epileptic 
patients without crisis. These results suggest that the epileptic crisis can be 
interpreted as a collective dynamical state characterized by the presence of a 
high degree of synchronization between different zones of the brain. 

A relevant finding of the present work is the presence of precursory maxima 
in both parameters r(t) and s(t) just before the onset of the epileptic crisis 
shown in Figure 4, while the corresponding EEG signal does not show appar- 
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ent changes at that moment. Our results suggest that the use of nonlinear 
dynamics techniques as tools to characterize EEG signals can contribute to a 
better understanding and diagnostic of some brain pathologies and, as we have 
shown, they may help in the prediction of epileptic events. Nonlinear dynam- 
ics techniques can also contribute to characterize other types of physiological 
signals. 
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